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1. Introduction 

Supersymmetry (SUSY) is one of the best studied extensions [|| of the Standard Model 
(SM). SUSY gives us also the benefit of introducing potential new sources of CP violation 
Q. As the small amount of CP violation in the SM is not sufficient to explain the baryon 
asymmetry of the universe, it is necessary to study all implications of the complex SUSY 
parameters. The search for SUSY will be one of the main goals of all future colliders. An 
e+e~ linear collider will be an ideal machine for the determination of the underlying SUSY 
parameters 

In this talk we present the results of our studies |Q, |5|, ^ of the effects of complex SUSY 
parameters on the phenomenology of the scalar top quark and scalar bottom quark system. 
Analysing the properties of 3rd generation sfermions is particularly interesting, because of 
the effects of the large Yukawa couplings. The lighter sfermion mass eigenstates may be 
among the light SUSY particles and they could be investigated at e+e~ linear colliders 
0, p. Previous analyses of the decays of the 3rd generation sfermions ii^2, bi,2, ^^1,2 and Uj- 
in the Minimal Supersymmetric Standard Model (MSSM) with real parameters have been 
performed in Refs. pj-pl]]. 

In the MSSM the SUSY parameters Af, fi and Mj (i = 1, 2, 3) are in general complex, 
where ^/ is the trilinear scalar coupling parameter of the sfermion fi, fi is the Higgs- 
higgsino mass parameter and Mi, M2 and M3 are the U(l), SU(2) and SU(3) gaugino 
mass parameters, respectively. We will first study the phase dependence of the decay 
branching ratios of the top squarks and bottom squarks, which are CP-even observables. 
As we will show, the phase dependence of these decay branching is indeed suitable to 
obtain informations about the SUSY CP phases. The situation is quite similar to that 
of the third generation slepton system, where the decay branching ratios of the staus fi^2 
and T-sneutrino Dr can be used to get information on the phases of the stau and gaugino- 



higgsino sectors |12|. In our study we will use the MSSM as a general framework and 
assume that the parameters At, Ab, fi and Mi have the phases ipAt, VAt, V^t and 9?u{i)i 
respectively (taking M2,3 real). We take into account explicit CP violation in the Higgs 



sector [^, |I4|, |T5|. Furthermore, we also take into account the constraints on the SUSY 
parameters which follow from the experimental data on the rare decay 6 — > S7 [^]. 

We will also consider a CP-odd observable in sfermion decays, which provides a more 
direct signal for the presence of CP phases. This observable is a CP-sensitive asymmetry 
which follows from triple product correlations ||6|, [l7[ |. 



2. Decay Branching Ratios of Top Squarks and Bottom Squarks 

Considering first top squark and bottom squark production, the reaction e~^e~ — > Qigj, 
Qi = ii,bi, proceeds via 7 and Z exchange in the s-channel. The tree- level cross sections ||8| 
of these reactions do not explicitly depend on the phases 99^ and ^pAq, because the Zqiiji 
couplings are real and in e~^e~ — > qiq2 only Z exchange contributes. The cross sections 
depend only on the mass eigenvalues 2 and on the mixing angles cos"^ 6q. Therefore, 
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they depend only implicitly on the phases via the cos((/?^ + (fAq) dependence of ruq-^ ^ 

In the following we will present numerical results for the phase dependences of the ti 
and bi partial decay widths and branching ratios. We will treat the fermionic decays 

iii^q' + Xk,i=h2,k = 1,2 (2.1) 



and 



and the bosonic decays 



and 



q^^q + xl^ = l,'^,k = l,..,A (2.2) 

q, ^ W^ + q'j,i = 1,2, j = 1, 2, j < i, (2.3) 

Qi^ H^ + q'j,i = 1,2, j = 1,2, j < i, (2.4) 

^2^^ + 91, (2.5) 

q2^ H, + qi,i = 1,2,3. (2.6) 



These partial decay widths depend on the SUSY parameters of the squark system Mq, 
Mjj, Mj^, tan/3, \^\, 93^, \At\, (pAt, \-^b\, VA,,, which determine the mass eigenvalues and 
mixing angles of the top and bottom squarks. In addition, in the chargino sector the SU (2) 
gaugino mass parameter M2 enters. The mass eigenvalues and mixing of the neutralino 
depend also on the U{1) gaugino mass parameter Mi with phase 93u(i) which, therefore, 
influences the partial decay widths of the top and bottom squarks into neutralinos. The 
main parameter of the Higgs sector is the charged Higgs boson mass mjj± , in addition to 
the parameters already introduced. 

We calculate the partial decay widths in Born approximation. In some cases the one- 
loop SUSY QCD corrections are important. The analyses of |ll|, |l^ suggest that a 
significant part of the one-loop SUSY QCD corrections to the partial widths of ti and 
bi decays (where the bottom Yukawa coupling is involved) can be incorporated by using 
an appropriately corrected bottom quark mass. In this spirit we calculate the tree-level 
widths of the ti and bi decays by using on-shell masses for the kinematic terms (such as 
phase space factors) and by taking running t and b quark masses for the Yukawa couplings. 
For definiteness we take mj"'^(mz) = 150 GeV, mp""^^" = 175 GeV, m^"'^(mz) = 3 GeV 
and TTiO^-^ii*^!' = 5 GeV. This approach leads to an "improved" Born approximation which 
takes into account an essential part of the one-loop SUSY QGD corrections to the ti and 
bi partial decay widths and predicts their phase dependences more accurately than the 
"naive" tree-level calculation. In the calculation of the CP violating effects in the neutral 
Higgs sector we take the program FeynHiggs-2.0.2 of [|l^]. 

In the numerical analysis we impose as theoretical constraint the approximate necessary 



condition for the tree-level vacuum stability |19]. Furthermore, as experimental constraints 



we take into account the mass bounds from LEP and Ap{t — b) < 0.0012 [21], as well 



as 2.0 X 10 ^ < B{b — > 57) < 4.5 x 10 ^ [16| assuming the Kobayashi-Maskawa mixing 
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Figure 1: Contours oiB{ti x\t) for tan/3 = 6, Ma = 300 GeV, = 500 GeV, V3u(i) = ^A, = 0, 
= 350 GeV, m^^ = 700 GeV, m^^ = 170 GeV, m^i = 600 GeV, with (a) \At\ = \Ab\ = 800 GeV 
and (b) ip^ — \Ab\ = \At\, assuming Mq > Mq. 



also for the squark sector. For the calculation of the b ^ width we use the formula of 
|22] including the 0(0^) corrections as given in |23|. We also take \Mi\ = 5/3tan^ 6'ly-^2• 
As a first example we show in Fig. |l] (a) the contour plot for B{ii Xit) ^ ^ function 
of if At and iff, for (m^-^ , m^-^ , J = (350,700, 170) GeV, tan/3 = 6, Ma = 300 GeV, \fi\ = 
500 GeV, l^tl = l^fel = 800 GeV, V7u(i) = fA^ = and mH± = 600 GeV, assuming Mq > 
Mjj. For the parameters chosen the ^pAt dependence is stronger than the (^^ dependence. 
The reason is that these phase dependences are caused mainly by the t^ -tpt mixing term, 
where the dependence is suppressed by cot /3. The dependence is somewhat more 
pronounced for ipAt ~ tt than for ipAt ~ 0, 2-ir. In Fig. |l| (b) we show the contour plot of 
Xit) as a function of ipAt and \At\ for ip^ = and \At\ = \Ah\. Clearly, the (pAt 
dependence is stronger for larger values of \At\. 

For the heavier top squark ^2 niore decay channels are open. In Fig. ^ (a) we show the 
branching ratios for t2 Xi,2^ and ^2 — > X2,3,4* ^ function of tpAt for (m^^ , m^-^ , m^^ ) = 
(350, 800, 170) GeV, tan /5 = 6, M2 = 300 GeV, = 500 GeV, | A| = \ At\ = 500 GeV,(p^ = 
(^U(^) = cp^^ = 0, m^^ = 350 GeV, m^^ = 800 GeV, m^^ = 170 GeV and m^^± = 350 GeV, 
assuming Mq > Mq. The ipAt dependence of B{i2 Xi'2^) is due to a direct phase 
effect, which explains that the shape of -B(t2 Xi2^) ^^^^ =t cos 99At)- Also the phase 
dependence of the branching ratios into neutralinos is mainly due to a direct phase effect. 
B{i2 — > ^ '^^^y weak phase dependence like (10 + cos ipAt)- In ^{^2 Xs^) the 

mixing phase enters, resulting in a shape like (1+cos (pAt) for the branching ratio. Similarly, 
Bih behaves hke (1 — cos ipAt ) ■ 

In Fig. ^ (b) we show the branching ratios for the bosonic decays i2 — ^ Zii and 
t2 — > Hiti [i = 1, 2, 3) for the same parameter values as above. The shape of B{t2 Zti) 
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Figure 2: ipA^ dependence of branching ratios of the decays (a) t2 (solid, black/gray), 

^2 — > X2/3/4* (dashed, black/gray/Ught gray) and (b) ^2 Zii (dashdotdotted), i2 ^^1/2/3*1 
(long dashed, black/gray/Hght gray) for tan/3 = 6, M2 = 300 GeV, = 500 GeV, \Ab\ ^ \At \ ^ 
500 GeV, ip^ = (/3u(i) = VA, = 0, m^^ = 350 GeV, = 800 GeV, mg^ = 170 GeV and 
mH± = 350 GeV, assuming Mq > Mjy. Only the decay modes with 5^1% are shown. The 
shaded areas mark the region excluded by the experimental limit B{b — > sj) < 4.5 x 10~^. 



is like (1 — cos ip At), which is caused by the 6^ dependence of the corresponding coupling. 
Quite generally, the phase dependence of r{t2 — > H^ti) is the result of a complicated 
interplay among the phase dependences of the neutral Higgs boson masses, the top squark 
mixing matrix elements, the neutral Higgs mixing matrix elements and the direct top 
squark-Higgs couplings. In the present example the fAt dependence of the partial widths 
r(t2 —>■ -ffi, 2,3^1) is mainly due to the ipAt dependence of the top squark mixing matrix 
and the squark-Higgs couplings, whereas the ipAt dependence of the neutral Higgs mixing 
matrix is less pronounced in this case. 

Coming now to the discussion of the decays of the bottom squarks 61^2, we show in Fig. ^ 
the partial decay widths and the branching ratios of 61 Xi 2^' W~ii as a function 

of LP At, for m^^ = 350 GeV, m^^ = 700 GeV, m^-^ = 170 GeV, tan/? = 30, mH± = 150 GeV, 
M2 = 200 GeV, = 300 GeV, \Ab\ = \At\ = 600 GeV, (/^^ = tt and (^a* = = 0, 

assuming Mq > M^. In the region 0.57r < (pA^, < I-Stt the decay 61 — > H~ti dominates. 
The ipAi, dependence of r(6i H~ti) is due to the behaviour of the squark-Higgs coupling. 
The partial decay widths r(6i — > Xi2^) ^^^^ almost ipA^, independent because the tpA^ 
dependence of the bottom squark mixing matrix nearly vanishes for tan (3 = 30. Hence the 
ipA^ dependence of the branching ratios B{bi Xi 2^) i^ caused by that of the total decay 
width. r(6i — > W~ti) is suppressed because hi ~ and ii ~ in in this scenario. 

We have also estimated what accuracy can be expected in the determination of the 
underlying MSSM parameters by a global fit of the observables (masses, branching ratios 
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Figure 3: ipAb dependences of (a) partial widths and (b) branching ratios of the decays bi — > }(^h 
(solid), 61 (dashed), 61 H~ii (dotted) and 61 —f W~ii (dashdotted) for tan/3 = 30, 

M2 = 200 GeV, = 300 GeV, \At\ = \At\ = 600 GeV, ip,, = tt, ipA, = ipv{i) = 0, m^^ = 350 GeV, 
m^^ = 700 GeV, m^^ = 170 GeV and toh± = 150 GeV, assuming Mg > M^,. 



and production cross sections) measured at typical linear collider experiments with po- 
larized beams. We have found that under favourable conditions the fundamental MSSM 
parameters except Af^ can be determined with errors of 1% to 2%, assuming an inte- 
grated luminosity of 1 ab~^. The parameter At can be determined within an error of 2 - 
3 % whereas the error of A^ is likely to be of the order of 50 %. More details can be found 
in @. 

3. CP Asymmetries in Sfermion Decays 

As the sfermions are scalar particles and they usually have two-body decay modes, it is 
not straightforward to construct a CP sensitive asymmetry involving their decays. In the 
case of top squarks one can use a three-body decay for the other sfermions one has 
to consider cascade decays. In our paper P we have considered the decay chain 

/ - /Xi ^fXiZ^ fx\ i I Uxi q q), (3.1) 

where i = e, fi, r, and q denotes a quark. We have defined a T-odd correlation for the 
leptonic decay 

Ofdd = P/-(P^ XP^-), (3.2) 

and for the hadronic decays as 

Oo'dd = P/-(P9XP5)> (3.3) 
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where p denotes the three-momentum of the corresponding fermion. We define the corre- 
sponding T-odd asymmetries as 



^ r(o^l>o) + r(o^l<o) 

which by CPT are also CP asymmetries. This CP asymmetry is similar to that proposed 



in 1 24], however, we calculate the asymmetry in the full phase space of the decay chain. For 
the measurement of or A\ it is necessary to be able to distinguish between the charges 
of (.'^ and £~ or q and q. In the case i = e, ^, r this should be possible experimentally on 
an event by event basis at an e'^e~ linear collider A'i^ will be measureable in the case 
of g = c, 6, where flavour reconstruction is possible [ p5[ . 

As an example we have calculated A^ for fi decay, considering the decay chain fi — > 
X2 ^ Xii Z ^ I I, ioT I = e,fj.,T. As input parameters we have chosen ruf^ = 
300 GeV, ruf^ = 800 GeV, tan/3 = 10, \Ar\ = 1000 GeV, ^pA^ = 0, = 800 GeV, 
M2 = 280 GeV, using the GUT relation \Mi\ = 5/3 tan^ ewM2 

In Fig. 1^ we show the contour lines for the branching ratio BR(fi t Xi ^ ^) = 
BR(fi ^ T X2) X BR(X2 ^ ^Xi) X BR(Z ^ ^ I) in the ipM^-V^^ plane for M2 = 280 GeV 
and \fi\ = 300 GeV. For BR(fi ^ r ^ ^) we always sum over £ = e,fi,T. We choose 

> since in this case the Ti-r-X2 coupling |o]^2l is larger, which implies a larger 
branching ratio BR(fi — > r X2) than for < Mj^. In a large region of the parameter 
space we have BR(x2 -^Xi) = 1) and we take BR{Z £ £) = 0.1. In Fig. ^ we show 
the ipMi and c/?^ dependence of A^. The value of A^ depends stronger on ipMi, than on 
iPfj_ . The sign of A^ is essentially determined by the sign of (pMi ■ 

The asymmetries A'^^^ can be calculated from A^ by replacing the leptonic Z coupling 
by the Zqq coupling. This gives 

^fj") ^ 6.3(4.5) X A{. (3.5) 

Based on our results for the asymmetry A^ in fi TX2 ~^ ^'ri^£~ and the branching 
ratio we give a theoretical estimate of the number of produced fi's necessary to observe the 
T-odd asymmetry. As an example we take the point denoted by • in Fig. ^, with = 7r/2 
and LpMi = ''r/2. For this point BR « 2.5 x 10~^ and \A^\ « 3 x 10~^. For the decay 
fi bbx^T, on the other hand, BR ~ 3.6 x 10^^ and |„4^| ~ 1.9 x 10^^. In this example the 
asymmeties A^'^ should be measurable at an e^e~ linear collider with ^/s = 800 GeV and 
an integrated luminosity of 500 for = 300 GeV. It is clear that detailed Monte 
Carlo studies taking into account background and detector simulations are necessary to get 
a more precise prediction of the expected accuracy. For a Monte Carlo study on a T-odd 



observable in neutralino production and decay see [26|. 
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Figure 4: Contour lines of the branching ratio for fi — > x\tU and asymmetry in the ipMi-'Pi_i 
plane for |^| = 300 GeV and Ma = 280 GeV, taking tan/3 = 10, Ar = 1000 GeV, rrif^ = 300 GeV, 
TOfa — 800 GeV for Af^ > M^^. The point denoted by • is for the theoretical estimate of the 
necessary number of produced fi's (see text). 
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